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Abstract 

In charged colloidal suspensions, the competition between square-well attraction and long-range 
Yukawa repulsion leads to various stable domains and Wigner supercrystals. Using a continuum 
model and symmetry arguments, a phase diagram of spheres, cylinders, and lamellae is obtained 
as a function of two control parameters, the volume fraction and the ratio between the surface 
tension and repulsion. Above a critical value of the ratio, the microphase cannot be supported 
by the Yukawa repulsion and macroscopic phase separation occurs. This finding quantitatively 
explains the lack of pattern formation in simple liquids because of the small hard sphere diameter 
in comparison with the size of macromolecules. The phase diagram also predicts microphase 
separation at zero value of the ratio, suggesting the possibility of self-assembly in repulsive systems. 
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Microphase separation is ubiquitous in soft matter systemsQ, Q, S, fl 3, B fl 3- 

For 

example, microphase separation in block copolymers results from mixing of two or more 
chemically different polymer segments ^ . Competition of hydrophobic and hydrophilic head- 
groups of amphiphiles leads to micelles in the water solution ^. The self-assembly processes 
in charged colloidal suspensions and protein solutions can lead to the formation of stable 
domains and Wigner supercrystals (see Fig. at low temperatures. The investigation 
of charged colloids and protein solutions is particularly interesting because they can be 
approximated as one-component systems with an effective isotropic pairwise interaction 
after averaging out the degrees of freedom of dispersing medium. This effective interaction 
is composed of a hard-core potential, a long-range electrostatic repulsion, and a short-range 
attraction. The screened Coulomb repulsion is described by a Yukawa potential 
Since more complicated potential forms qualitatively lead to the same phenomena 
for simplicity, we model the short-range attraction by a square-well potential. The overall 
pairwise interaction between two colloidal particles separated by r is given by 
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where a is the colloidal diameter. The attraction depth e and repulsion strength A are 
temperature dependent, and e is usually greater than the average thermal fluctuation. For 
convenience all the length variables in this letter are dimensionless in units of the screening 
length (. The pairwise and isotropic interaction in Eq. represents one of the simplest 
self-assembly systems that can be studied explicitly. 
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FIG. 1: Super structures of domains: A) bcc lattice of spherical domains, B) fee lattice of spherical 
domains, C) a cylindrical domain and the resulting 2D triangular lattice, D) a lamellar domain 
and the resulting ID lattice. 
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At low temperatures, shapes and sizes of colloidal domains are narrowly distributed, and 
their arrangements in space are highly ordered. A perfect superlattice composed of identical 
colloidal domains is an ideal reference (ground) state. In this letter we concentrate on domain 
patterns in the ground state. Entropic effects induced by thermal fluctuations will be studied 
elsewhere. Without thermal effects, we determine the most stable (optimal) shape, size, and 
superlattice (see Fig. ^ for ground-state domain patterns by globally minimizing the energy 
density. At low temperatures, this evaluation is simplified by a continuum approximation: 
Colloidal particles are closely packed inside a domain and the characteristic domain size 
is much larger than the particle size. Similar models have been used in the study of two- 
dimensional (2D) lipid domains^]. In our continuum model, the short-range attraction gives 
rise to bulk adhesion and surface energy. The density of the adhesive energy is a constant 
in the leading order for a given colloidal number density p or equivalently a given volume 
fraction, = vrpa^/G. The bulk adhesion does not affect the domain shapes and is not 
included in this letter. For a domain with area S, the surface energy is Us = 76", where the 
surface tension 7 is proportional to the attraction depth e in the lowest order approximation. 
Throughout this letter, 7 is assumed to be independent of domain shapes. The sum of the 
long-range Yukawa potential is separated into two parts. 
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where t/y^l'^) is the intra-domain repulsion for domain m, UY\''n,n) is the inter-domain 
repulsion between domains m and n, and is the total number of domains. Here pi is 
the colloidal number density within domains and larger than the overall number density p. 
The self energy of a single domain is the sum of the intra-domain repulsion and the surface 
energy. El = Us + Ui^^\ The sum of the inter-domain repulsions results in the lattice energy, 
U-L = (1/2) Uy\'^^ '^)- ground state, all the domains in one phase are identical 

so that the energy of each domain is given by i?tot = Ei + Ul- In our continuum model, 
the morphologies of domain patterns are thus determined by minimizing i^tot/'^d, where 
is the domain volume. 
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We now apply the continuum model to study various domain patterns shown in Fig. 
Spheres have the smallest surface energy and highest spatial symmetry. For a spherical 
domain with radius R, the intra-domain Yukawa contribution f/y ■* is simplified using the 
Fourier transform technique. The summation of the surface energy Us and the intra-domain 
repulsion results in the self energy of a spherical domain, 

El = 1^ [2R' + 3{a - 1)R^ + 3 - 3(1 + R)'e-'''] , (3) 
R-^ 

where Eq = Vd£o arises from a constant energy density = ^Pi^C^ ^iiid the domain volume 
fd = AnR^/S. As shown in Eq. (jS)), the competition between the surface tension and the 
Yukawa repulsion is described by a single control parameter, a = 'y/irplAC^. In the low 
density limit, domains are far apart from one another so that an individual domain can be 
treated as an isolated system and the lattice energy can be neglected. The self energy has a 
minimum at a finite radius for a < 1, whereas Ei is a monotonously decreasing function of 
R for a > 1. Thus ac = 1 is a critical point: Spherical domains with finite sizes are stable 
for a < Oc, whereas the phase separation occurs for a > ac- The self energy minimum Ei^^ 
and the associated radius Rra are plotted in Fig. |21 In the limit of strong repulsion (a ^ 0), 
we obtain two asymptotic forms, -Rm ~ v^l5a/2 and Ei „^/Eo ~ 3i^3a^/5. The size of 
the spherical domain grows as a increases. In the phase separation limit {a — > a^), the 
stable radius and energy minimum are asymptotically given by R^ ~ (1 — a/oc)"^^^ and 
Ei^m/ Eq ~ 2[1 — (1 — a/ac)^^'^ /y/S], respectively. These asymptotic relations can be tested 
experimentally. 

For a finite density (volume fraction), the balance of the self energy Ei and the lattice 
energy [/l from inter-domain repulsions leads to a 3D supercrystal, e.g., body centered cu- 
bic (bcc) and face centered cubic (fee) lattices of spheres (see Fig. [TJA. and B). The lattice 
energy [/l depends on the inter-domain repulsion and the lattice structure. Using the spher- 
ical harmonic expansion method, the inter-domain repulsion between two spheres separated 
by r is derived as 

U?\r) = ^ [(^ + +{R- l)e^ ' (4) 

The spatial periodicity of a Wigner lattice is related to the volume fraction. For example, 
the lengths of primitive cells for bcc and fee lattices (see Fig. ^A. and B) are given by 
a = r](p~^^^R, where we have r]{cc = (167r/3)^/^ for the fee lattice, and r^bcc = (87r/3)^/'^ for 
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FIG. 2: The optimal radius R^i and the energy minimum Ei^^ (-E-tot.m) of a spherical domain in 
the isolated case (dotted lines) and in fee lattices at finite volume fractions, where the dashed lines 
are for (j) = 0.1 and the solid lines are for (p = 0.3. 

the bcc lattice. Using the spatial periodicity, we sum the inter-domain repulsions of the fee 
(bee) lattice and obtain the lattice energy as 

Cl= ^[(fl+lK« + (ii-l)e«]^ 

Eexp{-r]^-^/^Rxfi/2) 
^- ■ 

The reduced distance in the above equation is given by xa = [2(^? + n| + n| + nin2 + 
^2^3, + n^ni)^^'^ for the fee lattice, and Xft = [3(n^ + + 7^3) — 2{nin2 + n2n^ + narii)]^/^ for 
the bcc lattice, where rti can be any integers except for that all the rti are zero. The total 
energy of a domain is the sum of the self energy in Eq. Q and the lattice energy in Eq. 
Minimization of the energy density for a given lattice using dfi{Etot/ Eq) = leads to the 
stable radius R^a- Figure |21 shows that both and the energy minimum -Etot,m increase 
with the volume fraction. By examining the subtle difference (~ 10~^) of -Etot.m between 
fee and bcc lattices, we obtain the phase diagram displayed in Fig. El It demonstrates that 
the fee lattice is usually more stable than the bcc lattice except for small values of a. Our 
results extend the previous studies on the fcc-bcc transition in Wigner particle solids with 
the Yukawa potential jisl. 

Dispersing medium surrounded by colloidal particles can be considered as cavities. For 
> 0.5, cavities are dispersed in the sea of colloidal particles so that various shapes and 
structures of cavities can form in the same way as domains. The intrinsic mirror symme- 
try between conjugate colloidal domains and cavities requires that the optimal shape and 
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FIG. 3: The phase diagram for fee and bee supererystals eomposed by spherieal domains. 

structure of cavities at 0(> 0.5) are the same as those of colloidal domains at 1 — 0. As a 
result, the cavity phase and the domain phase must be equivalent at = 0.5, which can be 
achieved by the lamellar shape and its topological variations. Following the same domain- 
cavity symmetry argument, lamellae are expected to be preferred in the phase separation 
limit (a a~). 

Next we investigate lamellar domains. For a lamellar domain with a finite width h, and 
infinite length and height {L oo), its surface energy is given by Us = 27(L^ + 2Lh) ^ 
2EQah~^, where Eq = v^eq = L'^heQ. Following Eq. Q, we obtain the intra-domain repulsive 
energy. 



In the ground state, all lamellar domains are parallel and form a one-dimensional supercrystal 
(see Fig. HP). To obtain the lattice energy, we calculate the repulsion energy between two 
parallel lamellar domains separated by r. 



At the volume fraction 0, the distance between two arbitrary lamellar domains is given by 
where n is a nonzero integer. The lattice energy for each lamellar domain is given 



[h-l + e~^] . 
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We determine the optimal width /im from minimization of Etot/Eo- Similar to spherical 
domains, lamellar domains with a finite width can exist only for a < Oc- In the limit of 
strong repulsion {a — 0) at a finite volume fraction, the lamellar width and the total energy 
are asymptotically given by /im ~ {/60a/(l — 0)^ and i?tot,m/-So = 20+ ^^9(1 — (jy^a"^ /{2(j)), 
respectively. By comparing -Etot.m of spheres and lamellae, we obtain their relative stability, 
as shown in Fig. |3] As expected, lamellar shapes are more stable than spheres when a 
approaches the critical value Oc, or when the volume fraction approaches 0.5. 

The spatial symmetry of spheres is the highest and that of lamellae is the lowest. It is 
natural to speculate that intermediate phases exist between these two limiting structures. 
One typical example is the 2D triangular lattice formed by cylindrical domains with the 
azimuthal symmetry (see Fig. ^C). For a cylindrical domain with a finite radius R and the 
infinite height {L oo), the surface energy is given by Us = 7(27ri?^ + 2ttRL) ^ 2EQaR~^ 
where Eq = v^^o = t^R^Leq. The intra-domain repulsion is derived from the Neumann 
addition theorem as 

=2Eo[l-2/i(i?)ifi(i?)], (9) 

where In{x) and Kn{x) are the modified Bessel functions of the first and second kinds, 
respectively. In this letter we consider the triangular lattice, although other lattice struc- 
tures or mixed structures are also possible We apply the Neumann addition theorem to 
calculate the effective repulsion between two domains separated by r as 

u'^\r) = ^E,[h{R)fKo{r). (10) 

The length of the primitive equilateral triangle shown in Fig. is related to the volume 
fraction as a = {2Tt / ^/^xpY^'^R. Using this relation, we calculate the lattice energy ?7l from 
the lattice summation of the inter-domain repulsion. By solving dptiEtot/ Eq) = 0, we obtain 
the optimal radius R^ and the energy minimum i^^tot,m- Similar to other shapes, cylindrical 
domains are stable for a < ac- 

Comparing -Etot,m calculated from different shapes and structures yields the global energy 
minimum and thus the optimal ground-state domain morphology. The central result of this 
letter is the phase diagram in Fig. IH which describes the shape transformation between 
spherical, cylindrical, and lamellar domains. The ratio a between the surface tension and 
repulsion, and the volume fraction (p, are the two control parameters. Finite size domains 
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FIG. 4: A phase diagram for spherical, cyhndrical, and lamehar shapes. The optimal structure 
formed by cavities for (p > 0.5 is the mirror of that by colloidal domains at 1 — 0. 

can be stabilized for a < ac = 1, whereas phase separation is observed for a > ac- Domain 
patterns stabilized in the small attraction limit {a — 0) suggest the possibility of self- 
assembly processes in repulsive svstemsjl^. The basic features of shape transformation 
are consistent with symmetry arguments. In the low density limit (0 0), spheres are 
preferred at small values of a, whereas lamellae are preferred as the system approaches phase 
separation. At larger volume fractions, structures in low dimensions (2D and ID) become 
increasingly more stable. Spheres are unstable for > 0.19, while cylinders are unstable 
for > 0.35. At the equal volume fraction of colloidal domains and cavities (0 = 0.5), 
only the lamellar phase is stable. The mirror symmetry between domains and cavities is 
used to produce the right half of the phase diagram for > 0.5. The cylindrical regime 
completely separates spherical and lamellar regimes, demonstrating that 3D spheres undergo 
a transformation to ID lamellae via 2D phases. Although we only compute cylinders, other 
intermediate shapes may exist. 

In this letter we predict the shape transformation (see Fig. of colloidal domain pat- 
terns in the ground state. Our study presents a simple and exactly solvable model system 
for understanding self-assembling phenomena based on a pairwise and isotropic potential. 
Ground-state domain patterns do not incorporate entropic effects induced by thermal fluc- 
tuations at finite temperatures. Temperature effects can be partially included in the current 
model by introducing the temperature-dependent surface tension 7(T) and domain density 
Pi(T). At higher temperatures, a more systematic treatment should involve the calculation 
of interphase free energies, where the distribution of shapes and distortion of structures are 
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accounted. Along this direction, the stabihty of domains with small distortions, the liquid- 
solid transition of particles within spherical clusters, and the formation of stable clusters at 
finite temperature are under investigation Q]. 

This work is supported by the NSF Career Award (Che-0093210) and the US Army 
through the Institute of Solidier Nano-technologies at MIT. 
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